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Abstract 

In the adiabatic and weak-modulation quantum pump, net electron flow is driven from one 
reservoir to the other by absorbing or emitting an energy quantum fuv from or to the reservoirs. 
In our approach, high-order dependence of the scattering matrix on the time is considered. Non- 
sinusoidal behavior of strong pumping is revealed. The relation between the pumped current and 
the ac driving amplitude varies from power of 2, 1 to 1/2 when stronger modulation is exerted. 
Open experimental observation can be interpreted by multi-energy-quantum-related processes. 
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I. INTRODUCTION 



Generally speaking, the transport of matter from low potential to high potential ex- 
cited by absorbing energy from the environment can be described as a pump process. 
The driving mechanics of classic pumps is straightforward and well understood-. The 
concept of a quantum pump is initiated several decades ago 2 - with its mechanism in- 
volving coherent tunneling and quantum interference. Research on quantum pump- 
ing has attracted heated interest since its experimental realization in an open quantum 

^^3,4,5,6,7,8,9,10,11,12,13,14,15 , 16.17.18.19.20.21.22.23 .24.25.26 .27.28.29.30 

The current and noise properties in various quantum pump structures and devices were 
investigated such as the magnetic-barrier-modulated two dimensional electron gas^, meso- 
scopic one-dimensional wire&>^, quantum-dot structures^^^ 1 *^, mesoscopic rings with 
Aharonov-Casher and Aharonov-Bohm effect^, magnetic tunnel junctions 1 -, chains of tunnel- 
coupled metallic islands^, the nanoscale helical wire22.,the Tomonaga-Luttinger liquid^ 6 -, and 
garphene-based devices^ 2 - 1 ^. Theory also predicts that charge can be pumped by oscillat- 
ing one parameter in particular quantum configurations 25 . A recent experiment 2 ^ based 
on two parallel quantized charge pumps offers a way forward to the potential application 
of quantum pumping in quantum information processing, the generation of single photons 
in pairs and bunches, neural networking, and the development of a quantum standard for 
electrical current. Correspondingly, theoretical techniques have been put forward for the 
treatment of the quantum pumps^ 2 ^ 2 ^ 7 -. One of the most prominent is the scattering 
approach proposed by Brouwer who presented a formula that relates the pumped current 
to the parametric derivatives of the scattering matrix of the system. Driven by adiabatic 
and weak modulation (the ac driving amplitude is small compared to the static potential), 
the pumped current was found to vary in a sinusoidal manner as a function of the phase 
difference between the two oscillating potentials. It increases linearly with the frequency in 
line with experimental finding. 

Although the quantum pump has been extensively discussed in literature, little attention 
was paid to experimentally observed deviation from the weak-pumping theory with only 
the first-order parametric derivative of the scattering matrix considered. We improved the 
scattering approach by expanding the scattering matrix to higher orders of the time and 
modulation amplitude, which enables us to go further in investigation of the problem. 



2 



II. THEORETICAL FORMULATION 



We start with the scattering matrix approach detailed by Moskalets et al— to describe the 
response of a mesoscopic phase-coherent sample to two slowly oscillating (with a frequency 
uj) external real parameters Xj(t) (gate potential, magnetic flux, etc.), 

X 3 (t) = X 0J + X^-e^-w) + X^e-^ l -^\ j = 1, 2. (1) 

Xqj and X w j measure the static magnitude and ac driving amplitude of the two parameters, 
respectively. The phase difference between the two drivers is defined as = tpi — y>%. The 
mesoscopic conductor is connected to two reservoirs at zero bias. The scattering matrix s 
being a function of parameters Xj(t) depends on time. 

It is assumed that the external parameter changes sufficiently slowly to validate an "in- 
stant scattering" description. To investigate the deviation from the small amplitude X w j 
limit, we expand the scattering matrix s(t) into Taylor series of Xj(t) to second order at 
Xqj with the terms linear and quadratic of X Ut j present in the expansion, 

s (t) « s (X 0>j ) + s- u e** + s +w e-^ + s 2 + s^e 2 ™ + s^e" 2 ^, (2) 

with 

' s±» = E XujeFwds/dXj, 
i=i,2 

< h= E Xlj&s/dXl (3) 
i=i,2 

s±2. = \ E Xl^dHjdX). 

\ 3=1,2 

It can be seen from the equations that higher orders of the Fourier spectra enter into the 
scattering matrix. As a result, both the nearest and next nearest sidebands are taken into 
account, which implies that a scattered electron can absorb or emit an energy quantum of 
Tujj or Ihio before it leaves the scattering region. In principle, third or higher orders in the 
Taylor series can be obtained accordingly. However, the higher-order parametric derivatives 
of the scatter matrix diminish dramatically and approximate zero. Numerical calculation 
demonstrates that even in relatively large amplitude modulation, their contribution is neg- 
ligible. 

The pumped current depends on the values of the scattering matrix within the energy 
interval of the order of max (fegT, 2hw) near the Fermi energy. In the low-temperature limit 
(T — ► 0), an energy interval of 2ftw is opened during the scattering process. 
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The mesoscopic scatterer is coupled to two reservoirs with the same temperatures T 
and electrochemical potentials fx. Electrons with the energy E entering the scatterer are 
described by the Fermi distribution function fo(E), which approximates a step function at 
a low temperature. Due to the interaction with an oscillating scatterer, an electron can 
absorb or emit energy quanta that changes the distribution function. A single transverse 
channel in one of the leads is considered. Applying the hypothesis of an instant scattering, 
the scattering matrix connecting the incoming and outgoing states can be written as 

k(t)=^2 SaP (t)a p (t). (4) 

P 

Here s a p is an element of the scattering matrix s; the time-dependent operator is a a (t) = 
J dEa a (E) e - lEt / h ) and the energy- dependent operator a a (E) annihilates particles with to- 
tal energy E incident from the a lead into the scatter and obey the following anticommutation 
relations 

[fit (E) , ap (E')} = 5 a(3 5 (E - E') . (5) 

Note that above expressions correspond to single- (transverse) channel leads and spinless 
electrons. For the case of many-channel leads each lead index (a, (3, etc.) includes a 
transverse channel index and any repeating lead index implies implicitly a summation over 
all the transverse channels in the lead. Similarly an electron spin can be taken into account. 
Using Eqs. (2) and (4) and after a Fourier transformation we obtain 

b a (E) = J2 [h,apap (E) + s 2 , a pap (E) + s- u>a pap (E + tuv) 

P (6) 
+s +u ,, a pa p (E — fko) + s_2w, a pap (E + 2fiw) + s +2w ^ap (E - 2tiw)] . 

The distribution function for electrons leaving the scatterer through the lead a is fa^ (E) = 
V i a {E)b a {E) s j , where (• • • ) means quantum- mechanical averaging. Substituting Eq. (6) 
we find 

f (out) {E) = £ + /q {E) + | . |2 /q {E + ?i£<;) 

\hu,ai3\ 2 fo (E-hu) + |s_ 2 ^,a/3| 2 fo (E + IDw) + \s +2w , a p\ 2 fo (E - 2Huj)] . 

The distribution function for outgoing carriers is a nonequilibrium distribution function 
generated by the nonstationary scatterer. The Fourier amplitudes of the scattering matrix 
|s- w , a /3| 2 (|s+aj,a/3| 2 ) is the probability for an electron entering the scatterer through the lead 
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(3 and leaving the scatterer through the lead a to emit (to absorb) an energy quantum 
fvjj and \s^2ui,ap\ 2 (1 3+2^,0,3 1 2 ) is that of the energy quantum 2hu process. |so,a/3 + ^2,a / 3| 2 is 
the probability for the same scattering without the change of an energy with the second- 
order term .§2,0/3 much smaller than the zero-order term so,a/3 in weak-modulation limit 
(X U j <C Xoj) and can be omitted therein. 

Using the distribution functions fo{E) for incoming electrons and f° ut (E) for outgoing 
electrons, the pumped current measured at lead a reads 

ip = 2^ r & {e) ha {e) ) - ^ {e) ^ {e) ^ dE - (s) 

Substituting Eqs. (7) and (3) we get 

= 2tt E X w j x X W j % ~gx~"gx~^ srn (fPh ~ Ph) 

,, (9) 
+f E XlA^f ^ sin [2 - ^ i2 )]. 

Quantum pumping properties beyond former theory based on first-order parametric deriva- 
tive of the scattering matrix are demonstrated in Eq. (9). By taking higher orders of the 
Fourier spectrum of the scattering matrix into consideration, double fuv energy quantum 
(or a 2fajj energy quantum) emission (absorption) processes coact with single fko quantum 
processes. In the weak-modulation limit, the second term in the right-hand side of Eq. (9) is 
small, which implies that double ftw quantum processes are weak and therefore not observ- 
able. As the ac driving amplitude is enlarged, this term increases markedly and contribution 
from double fko quantum processes takes effect. As a result, the dependence of the pumped 
current on the phase difference between two driving oscillations deviates from sinusoidal 
and changes from sin0 to sin 20, which is observed in experiment 3 . Moreover, the relation 
between the pumped current and the ac driving amplitude X u j is reshaped. It is also seen 
that the linear dependence of the pumped current on the oscillation frequency holds for 
multi-quanta-related processes. In the next section, numerical results of the pumped cur- 
rent in a two-oscillating-potential-barrier modulated nanowire are presented and comparison 
with experiment is given. 

III. NUMERICAL RESULTS AND INTERPRETATIONS 

We consider a nanowire modulated by two gate potential barriers with equal width L = 20 
A separated by a 2L = 40 A width well (see Fig. 1). The electrochemical potential of the 
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two reservoirs fi is set to be 60 meV according to the resonant level within the double- 
barrier structure. The two oscillating parameters in Eq. (1) correspond to the two ac driven 
potential gates Xi$ (t) — > (t) with all the other notations correspond accordingly. We 
set the static magnitude of the two gate potentials Z7o,i = Uq^ — Uq = 100 meV and the ac 
driving amplitude of the modulations equal U Wt i = = U u . 

In Fig. 2, the dependence of the pumped current on the phase difference between the 
two ac oscillations is presented. In weak-modulation regime (namely U u is small), sinusoidal 
behavior dominates. Here, three relatively large U u is selected to reveal the deviation from 
the sinusoidal dependence. (The magnitude of the pumped current mounts up in power- 
law relation as a function of U u as shown in Fig. 3. The sinusoidal curve for small U u 
would be flat and invisible in the same coordinate range.) It can be seen from the figure 
that the I p -(j) relation varies from sinusoidal (sin</>) to double-sinusoidal (sin 20) as the ac 
oscillation amplitude is increased. The interpretation follows from Eq. (9). The single fko 
quantum emission (absortion) processes feature a sinusoidal behavior while the 2fiw quantum 
emission (absortion) processes feature a double-sinusoidal behavior when the Fourier index is 
doubled. As U u is increased, double ftw quantum processes gradually parallel and outweigh 
the single ftw quantum ones. It is also demonstrated that when the single ftw quantum 
processes have the effect of sin0 dependence, the double fiw quantum processes induce a 
— sin 20 contribution with a sign flip, which can be understood from the sign change of the 
derivative of the scattering matrix. The effect of three- and higher luo quantum processes 
is small even for large comparable to Uq. The experimental observations 3 as a deviation 
from the weak-modulation limit are revealed by our theory. 

Experiment^ also discovered that for weak pumping the dependence of the pumped cur- 
rent on the pumping strength obeys a power of 2 relation, as expected from the simple 
loop-area argument^; for strong pumping, power of 1 and 1/2 relation is observed beyond 
former theory. We presented in Fig. 3 the numerical results based on our theory of the 
Ip-Uuj relation at a fixed <fi. To demonstrate its power-law dependence, natural logarithm of 
the variables is applied. From Eq. (9), it can be seen that for large ac driving amplitude U w , 
contribution of double Huj quantum processes (formulated in the second term on the right 
hand side of the equation) causes the I p -U w relation to deviate from its weak-modulation 
limit, the latter of which is I p oc U%. For different phase difference between the two ac 
drivers, the deviation is different. At <ft — n the pumped current is invariably zero regardless 
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of the order of approximation determined by time-reversal symmetry. At <fi = tt/2, sin 20 is 
exact zero, and no difference is incurred by introducing higher order effect. If we shift the 
value of <p to 0.497T, the abating effect of the double hio quantum processes has the order 
of with the small second-order parametric derivative of the scattering matrix smoothing 
that effect a bit. Consequently, a power of 2 — > 1 — > 1/2 relation is obtained and visual- 
ized by our curve fit, which is analogous to experimental findings. For different values of 
4>, sharper abating and augmental effect occurs with analogous mechanisms. It is possible 
that the experiment^ was done at the phase difference close to n/2 while trying to approach 
maximal pumped current in the adiabatic and weak-pumping limit. 

IV. CONCLUSIONS 

Based on the ac scattering approach, we go further to expand the time-dependent scat- 
tering matrix to higher orders of the modulation amplitude and the time. It is demonstrated 
in our theory that 2hu quantum emission (absorption) processes coact with those of single 
huj quantum when we go beyond the small-frequency and weak- modulation limit. Nonsinu- 
soidal dependence on the phase difference between two oscillating modulations is incurred by 
higher order Fourier components. The pumped current versus modulation amplitude rela- 
tion has a power law of 2 — > 1 — > 1/2 passage with the increase of the oscillating amplitude. 
Numerical results for a two-ac-gate modulated nanowire interpret experimental findings at 
large ac driving amplitudes. 
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FIG. 1: Schematics of the quantum pump: a nanowire modulated by two ac driven potential 
barriers. 
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FIG. 2: Pumped current as a function of the phase difference between the two modulations for 
different ac driving amplitudes. 
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FIG. 3: Pumped current as a function of the ac driving amplitude XJ W along with fits to I p oc U"^ (red 

1/2 

solid circle) below 35 meV, I p oc U u (green upward triangle) below 41 meV, and I p oc Uj above 
41 meV (blue downward triangle). To demonstrate its power-law dependence, natural logarithm 
of the variables is applied. The phase difference between the two ac driver <p = 0.497T. Inset is the 
zoom-in of the circled region. 
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